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Abstract 

Hill-determinant method is described and shown applicable within the so called PT- 
symmetric quantum mechanics. We demonstrate that in a way paralleling its tradi- 
tional Hermitian applications and proofs the method guarantees the necessary asymp- 
totic decrease of wave functions ip(x) as resulting from a fine-tuned mutual cancel- 
lation of their asymptotically growing exponential components. Technically, the rig- 
orous proof is needed/offered that in a quasi- variational spirit the method allows us 
to work, in its numerical implementations, with a sequence of truncated forms of the 
rigorous Hill-determinant power series for the normalizable bound states ip(x). 

PACS 03.65. Ge, 03.65.Fd 

1 Introduction 



One of the key sources of an enormous popularity of one-dimensional anharmonic- 
oscillator Schrodinger equations 

d 2 \ 

+V(x)\ ip(x) = Eif>(x), xe (-00,00) (1) 

with polynomial potentials, say, of the quartic form 

V(x) = ax 4 + bx 3 + cx 2 + dx, a = l (2) 

lies, definitely, in their methodically illuminating role of apparently "next-to-solvable" 
examples in quantum theory. 
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Thirty five years ago one of the characteristic results in this direction has been 
obtained by Bender and Wu [1] who made a highly nontrivial observation that in some 
models of such a type, all the pertaining bound-state energies E n may be understood 
as values of a single analytic function when evaluated on its different Riemann sheets. 
This observation (made at b = d = 0) was a real breakthrough in the intensive con- 
temporary perturbative analysis of eq. (1) + (2). Later on, it has been complemented 
by the alternative semi-classical [2] and quasi- variational [3] sophisticated treatments 
of eq. (2). In this sense, the model provided an extremely useful insight in the possible 
structures encountered within relativistic quantum field theory [4]. 

The hypotheses of ref. [1] have also been explicitly verified in the phenomeno- 
logically equally important cubic-oscillator limit a — > [5]. In parallel, the closely 
interrelated semi-classical and quasi- variational (often called Hill-determinant, HD) 
constructions of bound states were subsequently generalized to virtually all polyno- 
mials V(x). In the former semi-classical setting, an updated sample of references is 
offered by these proceedings. In contrast, the parallel HD approach to polynomial 
V(x) [6] seems "forgotten" at present, in spite of its extremely useful capability of 
complementing the perturbative and semi-classical constructions. This is one of the 
reasons why the forthcoming text has been written for the same proceedings. 

A nice illustration of the HD - semi-classical complementarity of the corresponding 
algorithms has been thoroughly described in ref. [7] (here, we shall skip this illus- 
tration referring to the original paper containing a few numerical tables and further 
references) . 

The HD method emphasizes certain important non-perturbative features of the en- 
ergies. Again we would like to recollect an elementary illustrative two-dimensional ex- 
ample of ref. [8] , or a more extensive one-dimensional sextic-oscillator study [9] where 
the HD-type energies have been specified as roots of a convergent continued fraction. 
Unfortunately, in the contemporary literature at least, the acceptance and credibility 
of the HD approach to polynomial V(x) has been marred, in the quasi- variational 
context as well as in perturbation theory, by a few unfortunate misunderstandings 
concerning the related mathematics. 

For a sketchy explanation of the latter point the reader should consult the twenty 
years old Hautot's concise review [10]. A more thorough account of the problem 
may be found in my own long and unpublished dissertation [11] including references 
to the misleading papers (responsible for the unfortunate continuation of the whole 
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misunderstanding up to these days) as well as a detailed account and an exhaustive 
specification of the domain of validity of the HD techniques in Hermitian cases with 
polynomial potentials. 

The present updated HD review is going to move one step further. The reason is 
that in the light of some most recent semi-classical analyses of anharmonic oscillators, 
another important motivation for a renewal of interest in HD philosophy might be 
found in the formal parallels between these two techniques, both having to deal with 
complexified coordinates and both being based on the infinite-series expansions of the 
energies E and of the related wave functions ip(x), respectively. It is also worth a 
marginal note that in both of them a different role is played by certain exponentially 
suppressed small contributions. 

2 VT— symmetric quantum mechanics 

Within the framework of the pure quantum physics, the idea of a complexification of 
the couplings or coordinates is in fact neither too popular nor formally welcome. A 
characteristic comment has been made by the author of [4] whose historical remarks 
clarified the situation in the seventies when people imagined that the convergence 
of perturbation series would imply that the "observables" may be well defined even 
when the Hamiltonians themselves become manifestly non- Hermitian. 

For this reason, perhaps, Carl Bender needed almost thirty years more to return to 
the complexified models in the two pioneering papers [12] where the two new versions 
of eq. (2) have been studied with the purely imaginary couplings b — i/3 and d = i5. 
The emphasis has been laid upon quite amazing an observation that the spectrum 
remains real in spite of the obvious fact that H ^ W . In the light of the possible 
applications of such an approach in Quantum Mechanics, the authors characterized the 
underlying specific non-Hermiticity as a "weakened Hermiticity" or '"PT symmetry" 
of the Hamiltonian. Indeed, in their most elementary models the operators V and T 
denoted the parity and the time reversal (i.e., in effect, the Hermitian conjugation), 
respectively. 

The methods they used were mainly perturbative and semi-classical. In what 
follows, we intend to answer the question whether some models of VT — symmetric 
quantum mechanics could admit the constructive solution of their Schrodinger equa- 
tions by the techniques of the HD type. For the sake of definiteness and in a way 
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employing our recent results in [13], we shall mainly review the results concerning the 
above-mentioned quartic problem 



again, with the real values of /3, c and 5 and, therefore, with the manifestly non- 
Hermitian, VT— symmetric form of the Hamiltonian. 

Before getting in more details, let us only add that during the rapid development 
of the field (summarized, e.g., in the proceedings of the recent Workshop which have 
to appear in the October issue of Czechoslovak Journal of Physics this year), the 
"parity" V should be understood, in a broader sense, as an arbitrary (Hermitian and 
invertible) operator of a pseudo-metric rj/ pseu( io) in Hilbert space [14]. Similarly, the 
symbol T should be understood as an antilinear operator mediating the Hermitian 
conjugation of H [15]. 

In the context of physical applications, the concept of the "parity" itself becomes 
insufficient and must be complemented by a "quasi-parity" operator Q [16] which 
makes the scalar products (between eigenstates of H ^ W) positively definite [17]. 
This means that we must perform a highly nontrivial transition from the non-singular 
and, by assumption, indeterminate original pseudo-metric V = r)( pseuc io) = vlpseudo) t° 
some other, positively definite "true" metric QV = r] + > 0, the knowledge of which 
allows us to call our Hamiltonians H ^ W "quasi-Hermitian" and "physical" again 
(cf. the reviews [18] for more details). In the other words, only the change of the 
metric V — > QV = r]^ quasi ) = m quasi , (where, at present, the name "charge" for Q = C 
is being preferred [19]) makes the non-Hermitian VT— symmetric models with real 
spectrum fully compatible with the probabilistic interpretation and other postulates 
of quantum mechanics in a way illustrated very recently on the non-Hermitian square 
well [20, 21]. 

3 HD construction 

The wave functions ip(x) defined by eq. (3) are analytic at all the complex x G C. In 
the asymptotic region of the large \x\ 3> 1 the semi-classical analysis reveals that 



dx 2 



+ x 4 + i(3x 3 + cx 2 + i5x) i/j(x) = Etp(x), 




(3) 




c p hy S e- x3/3 + c unphys e +x3/3 , < | Im (x) | < + Re {x)/y/Z, 
d P hy S e +xS/3 + d unphys e- x3/3 , < | Im (x) \ < - Re {x)/y/Z. 



(4) 
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The HD method proceeds, traditionally, in an opposite direction and employs the 
power-series formulae for ip(x) near the origin. Thus, in the VT— symmetric manner 
we write 

oo 

^ ansatz \x) = e- sx2 Y, h n (™) n , x e {-oo, oo) (5) 



n=0 



with a suitable (optional) value of s. The key reason may be seen in the success 
of such a strategy in the harmonic-oscillator case where such an ansatz converts the 
differential Schrodinger equation (3) into solvable recurrences. Here, of course, the 
resulting recurrent rule 



A n h n+2 + C n h n + 5 h n -i + 6 h n _ 2 — (3 /i n _ 3 + /i„_ 4 — 
A n = (n + l)(n + 2), C n = Asn + 2s-E, 6 = As 2 - c 



(6) 



is not solvable in closed form. At all the parameters (including also the unconstrained, 
variable energy E), it only defines the coefficients h n as superpositions 



h n — ho<7 n + h\LJ n , 



(7q = tO\ — 1, <7i = UJ = 



(7) 



where all the three sequences h n , a n and oo n satisfy the same recurrences with different 
initial conditions. Still, we may write the latter two funtions of n in a very compact 
determinant al form 



(Jn+i = (-iy 



det £„_ 



n\{n + l)\ ' 
with (m + 1)— dimensional matrices 



_ ]n det 

n\(n + 1)! 



n=l,2,... 



fC A 

5 A 1 

C 2 A 2 

1 ! C 3 ••• 

-p ■•. ■•. 

1 i ••. s 

V '•• ••• 



o 

C m -i 
5 



A m -2 





A m -1 

i 



(8) 
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(note: the second column of the matrix of the linear system (6) is omitted here since 
<j\ = 0) and 



a 



/ 





A 






\ 




Ci 





A l 








5 


c 2 





A 2 








5 


c 3 








1 




5 





A m -2 










Cm-l 


A m _! 


V 






1 


5 


c m o J 



(9) 



(now the first column is omitted in the light of our choice of c^o = 0). 

It remains for us to specify the norm p = \Jh,Q + hf, the ratio hi/h = tan£ and 
the physical energy E which enters the (m+ 1)— dimensional matrices (8) and (9). In 
the other words, our wave functions ip(x) must be made "physical" via the standard 
boundary conditions 



^ ansatz \X R ) = = ^ ansatz \-X L ), X L >1. 



(10) 



These conditions play, obviously, the role of an implicit definition of the above free 
parameters. Nevertheless, our ambitions are higher and in the HD approach people 
usually try to replace eq. (10) by a quasi-variational finite-matrix truncation of eq. 
(6) [10], i.e., equivalently, by innovated difference-equation boundary conditions 



h N = h N+1 = , iV>l. 



(11) 



with a much simpler numerical implementation. Of course, such a replacement of eq. 
(10) by eq. (11) is not always well founded [10] so that in each particular construction 
its validity must be based on a rigorous mathematical proof, in a way illustrated in 
what follows. 



4 The proof 

4.1 The s— dependence of the asymptotics of h n 

In the first step of the proof, the recurrences (6) should be read as a linear difference 
equation of the sixth order. The sextuplet of its independent asymptotic solutions 
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h n may be found by the standard techniques. Thus, the leading-order solution may 
be extracted from eq. (6) by its reduction to its two-term dominant form of relation 
between h n+2 and /i n _4- Thus, we put 

(31/3)n r(l+n/ 3) W 

where the integer p = 1, 2, . . . , 6 numbers the six independent solutions and where 
X(p) = exp[i(2p — l)7r/6]. The new coefficients g n = g n {p) vary more slowly with 
n. Thus, our first conclusion is that at the large indices n, all solutions decrease as 
h n ~ 0(n- n/3 ) at least. This means that in a way confirming our expectations the 
radius of convergence of our Taylor series (5) is always infinite. 

On this level of precision the size of our six independent solutions remains asymp- 
totically the same. In order to remove this degeneracy we amend equation (6) and 
having temporarily dropped the p's we get 

_ 4sA 4 (3\ ( g n \ 

9n+2 ~ 9n-4 ~ 9n ~ ^3 9n-3 + O J . (13) 

The smallness of 1/n 1 ^ 3 in the asymptotic region of n 3> 1 enables us to infer that 

9n = e ^ 2/3+ °(" 1/3 ), 7 = 7 (p) = s\ 4 (p) - (3 A(p)/4, (14) 

i.e., 

Re 7 (1) = Re 7(6) = ~/3 - |, Re 7 (2) = Re 7 (5) = s, 

Re 7 (3) = Re 7 (4) = ^-^. (15) 
We may conclude that whenever we satisfy the condition 

s > & (16) 

the general, six-parametric form of the Taylor coefficients 

K = Y,G p K{p) (17) 
p=l 

will be equivalent to the two-term formula in the leading order, 

h n = G 2 h n (2) + G 5 h n (5), n»l, (18) 
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i.e., only its two dominant components remain asymptotically relevant while, in the 
leading order, we may simply put G\ = G 3 = G 4 = G 6 = in eq. (17) at n >> 1. 
The same argument implies that for 

we get 

h ( G 3 h n (3) + G 4 h n (4), (3 > 0, n > 1, 
n \ GiMl) + G 6 M6), /3<0, n»l. 

while the degeneracy of more than two solutions survives at (3 = or at s = > 0. 

4.2 Single-term dominant exponentials in ijj( ansatz ) (x) 

We know that the shape of the function ^j( ansatz ) ( x ) is determined by the energy E 
and by the not yet fully specified choice of the two coefficients h and h\. In the 
language of the preceding section, each choice of the energy E and of the initial h 
and hi will generate a different, x— and n— independent pair of the coefficients G 2 
and G5 in h n at the sufficiently large s (the discussion of the smaller s will be omitted 
for the sake of brevity). 

As long as our aim is a completion of the proof of the validity of the replacement of 
the standard boundary conditions (10) (containing infinite sums) by the more natural 
HD truncation of recurrences due to eq. (11), we just have to parallel the Hermi- 
tian considerations of ref. [22]. Firstly, we remind the reader that at an arbitrary 
finite precision we always have E 7^ E (physical). This means that our infinite series 
^( ansatz )(x) as defined by equation (5) will always exhibit an exponential asymptotic 
growth in a certain wedge-shaped vicinity of the real line [12]. 

The main point of the HD proof is that due to the standard oscillation theorems the 
last (left as well as right) nodal zeros will move towards their (left and right) infinities 
with the decrease of E > E(physical) (say, at a fixed pair /io,i)- In the language 
of coordinates this means that the values of ip( ansatz ) (x) (both at some x > Xr and 
at another x < Xl) will change sign for E somewhere in between the (smaller) 
E (physical) and (larger) E(numerical). Such a simultaneous change of the sign of 
ip( ansatz ) ( x ) at any sufficiently large absolute value of the coordinate \x\ ^> 1 should 
be understood as resulting from our appropriate asymptotic estimate of ^j( ansatz ) ( x ) . 
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Although this sounds like a paradox, this estimate will always lead to an asymp- 
totic growth (since E ^ E physica i with probability one) so that it cannot be changed 
by the (key) omission of any finite number N of the exponentially small components 
~ exp[— sx 2 + 0(\nx)}. They may be safely ignored as irrelevant. This means that we 
may choose iV 3> 1 so that just the asymptotically dominant coefficients will play the 
role. Inserting (12) and (18) in ^( ansatz )( x ) ~ exp(-sx 2 ) En=N+i h n (ix) n we finally 
get 

G 2 X n (2)g n (2)+G 5 X n (5)g n (5) 



ansatz) 



x) ~ e 



E 

n=N+l 



{ix)\ 



\x\ > 1. 



(3 1 /3)nr(l + n/3) 

The validity of this formula is a strict consequence of the specific constraint (16) 
imposed (say, from now on) upon the admissible quasi- variational parameter s. 
Once we split ifj( ansatz ) ( x ) = ijj( ansatz \G 2 , G 5 , x) in its two components 

7 (2) n 2 / ? > + O^ 1 ' 3 ) 



^ ansatz \G 2 ,0,x) ~ G 2 e 



^ ansatz \0,G 5 ,x) ~ G 5 e 



E 1 

n=N+l 

oo 

n=N+l 



-X) 



cxp 



x exp 



(3 1 /3)«r(l + n/3) 
(5) n 2 ' 3 + Oin 1 ' 3 ) 



(3V3)n r ( 1 + n /3) 

we may apply the rule e z ~ (1 + z/t)*, t ^> 1 in the error term and get 



ij ( - ansatz \G 2 ,0, -y) 



and 



exp(— sy A 



tp( ansatz \0,G 5 ,y) 



OO J 

° 5 n S +1 (3V3)nr(l + n/3) 



1/ • 



1 + 



ATI/3 



1 + 



i 



ATl/3 



exp{—sy*) 

This is valid at all the large arguments y. Along the positive semi- axis y 1, 
both the right-hand-side summands are real and positive. They sum up to the same 
function exp[y 3 /3 + 0(y 2 )}. This is a consequence of the approximation of the sum 
by an integral and its subsequent evaluation by means of the saddle-point method. 
The same trick was used by Hautot, in similar context, for the V— symmetric and 
Hermitian anharmonic oscillators [10]. 

In contrast to the Hautot's resulting one-term estimates of ip, the present asym- 
metric, VT— invariant construction leads to the more general two-term asymptotic 
estimate 



(ansatz) 



(G 2 , G 5 , x) ~ G 2 exp[-a; 3 /3 + 0(x 2 )] + G 5 exp[x 3 /3 + 0(x% \x\ > 1. 
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As long as we deal with the holomorphic function of x, this estimate may be analyti- 
cally continued off the real axis of x. Near both the ends of the real line and within 
the asymptotic wedges |Im x|/|Re x\ < tan7r/6 we simply have the rules 

^ ansatz \G 2 ,G 5 ,x) ~ G 2 exp[-a; 3 /3 + 0{x 2 )}, Re x < -X L < -1 (21) 

and 

ij {ansatz \G 2 ,G 5 ,x) ~G 5 ex P [x 3 /3 + 0(x% l<^X R <Rex. (22) 

They are fully compatible with our a priori expectations and represent in fact the 
main step towards our forthcoming completion of the rigorous proof of the validity of 
the HD matrix truncation (11). 

4.3 The changes of sign of ^ anaatz ){x) at \x\ > 1 

Our complex differential Schrodinger equation (3) becomes asymptotically real, in the 
leading-order approximation at least. In a suitable normalization the wave functions 

^(ansatz) 

(x) may be made asymptotically real as well. Near infinity they will obey the 
standard Sturm Liouville oscillation theorems. As we explained, after a small decrease 
of the tentative energy parameter E > E(physical) the asymptotic nodal zero Xr or 
— Xl originating in one of our boundary conditions (10) will move towards infinity. 
This may be re-read as a doublet of conditions 

G 2 = G 2 (E , Co) = 0, G 5 = G 5 (E , Co) = (23) 

using an assumption that Co ~ ({physical) in the suitable parametrization of ho = 
p cosC and hi — p sinC with C £ (0, 2ir) at a convenient normalization p — 1. 

In the limit iV — > oo the two requirements (23) may be interpreted as equivalent 
to the truncation recipe (11). Indeed, at a fixed # >1 we may re-scale 

_ \ N (p) exp[ 7 (p)iV 2 / 3 ] 
h ( 3 i/3)Jv r(1 + iV/ 3) ' P 

where f p (E, Co) ~ F p ■ (E — E ). This enables us to write 

h N » (F 2 + F 5 )(E - E ) + 0[(E - E )% 
(N + 3y/ 3 h N+1 » [F 2 A(2) + F 5 \(5)](E - E ) + 0[(E - E Q f] 

due to equation (12). We see that the two functions G 2 , G 5 are connected with 
the Taylor coefficients h^ = h^{E Q , Co) an d frjv+i = h N+ i(E , Co) near the physical E 
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and Co by an easily invertible regular mapping. This means that the implicit algebraic 
boundary conditions (23) are strictly equivalent to the fully explicit requirements 

h N (E o ,{ ) = 0, h N+1 (E , Co) = 0, iV>l. (24) 

This completes our proof. 



5 Discussion 

5.1 Truncated secular Hill determinants 

By construction, our present result may be re-read as a demonstration that the intu- 
itive quasi-variational square-matrix truncation of our recurrences represents in fact 
a mathematically well founded approximation recipe. The evaluation of both the en- 
ergies and wave functions may be started at any approximative cut-off N < oo and 
the solution of the linear algebraic problem 



Co 





Ao 


5 


Ci 





e 


5 




-(3 


e 




1 


-(3 





\ 



A 1 



V 



1 -p 9 



An-3 



CW-i J 



I 



ho 
h 2 



\ 



h-N-3 
h-N-2 

V h N -! J 



(25) 



must only be completed by the limiting transition iV — > oo. As long as the energy 
only enters the main diagonal, C n = 4sn + 2s — E, we may determine all the ap- 
proximate low-lying spectrum by the routine N x iV— dimensional asymmetric-matrix 
diagonalization. With s = 2 and a = c = (3 = 5=l, Table 1 illustrates the practical 
implementation as well as the highly satisfactory rate of convergence of the numerical 
HD algorithm. 



5.2 Outlook 

We have re-confirmed that the HD methods bridge a gap between the brute-force 
variational algorithms and sophisticated analytic semi-classical constructions. In the 
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present context, we emphasized that the HD techniques may be combined not only 
with the complexification of the variables (characteristic for the latter approach) but 
also with the requirements of a numerical efficiency (pursued, usually, in the former 
setting). 

In the conclusion we may note that another link connects also the HD and pertur- 
bative techniques. In the latter, semi-analytic type of considerations, an indispensable 
role is played by asymptotic expansions, say, of the power-series form 

Z(X) = '£X k z k + R M (X) (26) 

k=0 

where A represents a - presumably, "sufficiently" small - coupling constant while the 
typical example of the observable Z(X) is a bound-state energy. In the present HD 
setting we employed a similar ansatz (26) where the variable A denoted the coordinate 
x while the function Z{X) represented the wave function. 

Currently, a numerical credibility of the similar perturbative (and also semi- 
classical) constructions is being enhanced via an improved treatment of the error 
term Rm(X). In general, this does not seem to be the case in the present HD context. 
The reason is that although for the bound states the exponentially growing right-hand 
sum must be identically zero, we are still not interested in the explicit evaluation of 
the exponentially suppressed remainder Rm{X) (which would be hopeless) but merely 
in the bracketing property of the wave-function sum J2h=o X k Zk 7^ evaluated slightly 
below and slightly above the actual physical energy. 

In this sense, we may summarize that the HD "trick" is different and enables us 
to pay a more detailed attention to the new perspectives opened by the consequent 
complexification of A. Here, such a perspective enabled us to extend the standard 
mathematical background of the Hill-determinant constructions of bound states. We 
have seen that the old principles stay at work, making use of the same mechanism 
as in the standard Hermitian cases, viz, the constructively employed asymptotic can- 
cellation of the growing exponentials in the exponentially decreasing wave functions 
ip(x). 
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Table 



Table 1. The convergence of HD energies. 





energies E n 


n 





1 


2 


3 


4 


N 












15 


1.69347 


5.106 


9.152 


13.043 


17.82 


20 


1.691638 


5.1256 


9.2800 


14.050 


19.24 


25 


1.691579 


5.12344 


9.2581 


13.869 


18.79 


30 


1.691590 


5.12361 


9.2617 


13.883 


18.89 


35 


1.691590 


5.12358 


9.2615 


13.879 


18.88 
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